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Abstract. The Kuramoto model when considered over the full space of phase angles 
[0, 27t) can have multiple stable fixed points which form basins of attraction in the 
solution space. In this paper we illustrate the fundamentally complex relationship 
between the network topology and the solution space by showing that determining the 
possibility of multiple stable fixed points from the network topology is NP-hard for 
the weighted Kuramoto Model. In the case of the unweighted model this problem is 
shown to be at least as difficult as a number partition problem, which we conjecture to 
be NP-hard. We conclude that it is unlikely that stable fixed points of the Kuramoto 
model can be characterized in terms of easily computable network invariants. 


PACS numbers: 05.45.Xt 

1. Introduction 

The Kuramoto model [16] was originally motivated by the phenomenon of collective 
synchronisation whereby a system of coupled oscillating nodes will sometimes lock on 
to a common frequency despite differences in the natural frequencies of the individual 
vertices. This model has since been applied to a wide variety of application fields 
including, biological, chemical, engineering, and social systems (see the survey articles 
including [1], [2], [22], and recently [8]). While Kuramoto studied the infinite complete 
network it is natural to consider finite networks of any given topology. This would 
correspond to a notion of coupling that is not universal across all node pairs, but rather 
applies to a subset of all possible links. For example the work patterns of human 
individuals in an organisation might enjoy a coupling effect in relation to pairs of 
individuals that have a working relationship. In this article we explain the significance 
of stable fixed point solutions to the model, and explore the algorithmic complexity of 
finding such solutions. 
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1.1. The model 

Each node has an associated phase angle 9 i: , as well as its own natural frequency a;*. 
The basic governing equation is the differential equation: 

n 

6i =uJi + k y Ajj sin (fl 3 - - 9i),i = 1, .., n. (1) 

1 = 1 

where A is the adjacency matrix of the network and k is a coupling constant which 
determines the strength of the coupling. Note that each 9i is understood to denote a 
function of time t. An important special case is where all the natural frequencies cuj 
are equal. We shall call this case homogeneous and the general case as inhomogeneous. 
Typically A is a (0,l)-matrix, however it may also be any real positive symmetric matrix 
reflecting coupling constants that vary for each edge. This case is of particular interest 
to us and we specify this by referring to it as the weighted Kuramoto model. It is well 
known that results for unique stable fixed points can be obtained for restricted phase 
angles 9i, say in the range [0,7r]. In this article we are concerned with results about 
stable fixed points over the full range [0,27 t). 


1.2. Synchronisation 


It has been observed that for arbitrary initial phases, networks synchronize in that some 
of the node phases converge to the same, or nearly the same phase angle, while the phase 
frequencies 9{ converge to a common value. Meanwhile, the remaining nodes behave 
non-uniformly or ’drift’. Moreover, at a sufficiently large critical coupling constant k, it 
becomes possible for all nodes to participate in synchronized behaviour. By summing 
over equation (1), the sine terms cancel and the average frequency of the nodes is a 
constant ui so that 


1 

n 




1 

n 



( 2 ) 


3 =1 3 = 1 

Define a frequency fixed point as a situation in which all the node frequencies are equal 
and are fixed over time. By equation (2), this is characterized by 


9i = u,i = 1, 


n. 


( 3 ) 


It follows that for a frequency fixed point all phase angles remain constant. We note 
that a frequency fixed point may also be accompanied by phase synchrony whereby the 
nodes have the same, or nearly the same, phase angles. This is however not necessarily 
the case, and the phase angles of a frequency fixed point may be significantly different. 
Examples of this kind are the ring networks [19]. Combining equations (1) and (3), the 
phase differences at a frequency fixed point satisfy 

n 

Co = LUi + k Aij sin(6 l J - — 9f),i = 1,.., 

3 = 1 


n. 


( 4 ) 
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Henceforth, we shall use the term fixed point to refer to a frequency fixed point satisfying 
equation (4). The character, number and location of the fixed points of a network are 
clearly important in understanding the types of dynamics that are possible from all 
possible initial phases of the system. In particular, understanding the nature of all 
stable fixed points is particularly important, since each represent an attractor set of 
positive n-dimensional volume within the set of all states. Fixed points can have a 
range of attractor types ranging from unstable single points, 777-dimensional ( 777 < n) 
saddle point structures that are partially stable, to 77 -dimensional volumes for stable 
fixed points. 

1.3. Stability 

Let {0*,i = 1, .., 77 } be a fixed point so that by equation (3) 

9*(t)=Q,i = l,..,n. (5) 

Since 9*(t)—9*(t ) is a constant for all i,j we use the compact notation 9*{t)—9*{t) = A* t . 
A linear stability analysis shows the fixed point to be stable (see for eg. [24]) if and only 
if the matrix M is negative definite where 

n 

M tj = A i:j cos(AL) - S i:j ^2 A ik cos (Afc;), (6) 

k =1 

and where the Kronecker delta function is defined as 8ij — 1,0 for i — j and i 7 ^ j, 
respectively. For a fixed point 9 *, i = 1, .., 77 , we can rearrange equation (4) as 

_ n 

'A—Al — ^2 A.ij sin(A*-), i = i, 77 . (7) 

7=1 

It follows that as the coupling constant increases in relation to the natural frequencies 
the left hand side of equation (7) tends to 0, in the limit giving: 

n 

^2 Aij sin(A7) = 0, i = 1, 77. (8) 

7=1 

Note that this limiting equation is the same as that obtained in the homogeneous case, 
and also that this equation is independent of the value of the equal natural frequency. 
Thus as the coupling constant increases the fixed point for unequal natural frequencies 
converge to the solutions for the homogeneous system of equations (8). Furthermore, in 
the homogeneous case it is sufficient to consider the case where the natural frequencies 
are all zero. For this reason the fixed point character for the case where all natural 
frequencies are zero is fundamental to understanding the limiting behaviour of the 
general inhomogeneous case. Similarly the stability condition does not depend on the 
natural frequencies (see equation (6)). Thus the homogeneous case can in some sense 
be considered as a guage of the contribution of the purely topological effects (from A^) 
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on fixed points and stability. Though the relationship between the network topology 
Aij and the stable fixed points is a complex one, we have the following useful necessary 
condition [24], Let {&*, i = 1, n} be any stable fixed point solution to the homogeneous 
Kuramoto model. Then for any proper node subset X we have the equality 

Ai i sin ( A P = °- ( 9 ) 

(i,j)e(x,x°) 


and inequality 

cos(A’ ; ) > 0. (10) 

p,j)e( x,x°) 

It follows that if for every edge (i,j) if |A*J > 7 t/ 2, then the fixed point &*, i — 1,.., n is 
unstable. It is also well known [19] that if for each link (i,j) A* ? < 7t/ 2 then the fixed 
point is stable. These equations are useful in resolving the mixed case where some |A*J 
are less than 7 t/ 2 and some are greater than 7 t/ 2. We shall often find the case where 
the sets X in (9) and (10) are single nodes particularly useful, where (9) and (10) take 
the form 

n 

^2 A ij sin(AL) = 0, i = 1,.., n. (11) 

3 = 1 
n 

Y A ij cos(AL) > 0, i = 1,.., n. (12) 

j =i 

It is useful to visualize the phase angles as points on the unit circle within the complex 
plane where the nodes are positioned around a circle, which we shall refer to as a phase 
angle diagram. Thus we let z^i = l,..,n be given by Zj = e l0 > = co s(9*) +isin(9*). 
This is illustrated in figure 1 with the representation of the hexagonal ring graph which 
forms a non-zero stable fixed point. There is also a simple physical intuitive way of 
understanding the phase angle diagram by interpreting each edge as an elastic link 
between beads on a frictionless ring. A fixed point corresponds to the elastic forces at 
each bead balancing out, and stability is manifested by any small perturbations of the 
beads returning to the fixed point (see [8]). 


2. Non-zero stable fixed points 

Though the nature and presence of stable fixed points is well understood for certain 
classes of graphs, general results for arbitrary graph topologies seem difficult to obtain. 
Thus complete graphs can have only one single stable fixed point [18] as can all dense 
graphs with degree at least 0.9395(n — 1) [24], Also trees can have only one stable 
fixed point [4] though the nature of all fixed points and the relationship to the coupling 
constant is well understood [7]. On the other hand ring graphs do have multiple stable 
fixed points the nature of which is also well understood [19] . In contrast we show that the 
most basic of questions about the stable fixed points of the weighted Kuramoto model 
is unlikely to be computable efficiently. Specifically, determining the presence of any 
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Figure 1 . Representation of nodes on the complex plane 


non-zero stable fixed points from the network topology A,j even for the homogeneous 
system is NP-hard. This means that an efficient (polynomial time) algorithm to solve 
this problem is either very difficult to find or does not exist (see the standard guide 
to algorithmic complexity [10]). This remains true even if the phase angle differences 
across edges ^ 0) are restricted to A*J < tt/2. For the unweighted case we show 
that determining the presence of any non-zero stable fixed points from the network 
topology Aij is at least as difficult (up to polynomial complexity) as solving a particular 
number problem. To establish the former we shall refer to the Partition problem [10]. 
In terms of related algorithmic complexity results we are aware only of the result of 
[20] which shows that checking positive definiteness of symmetric interval matrices is 
NP-hard. This problem is significantly different to the problem presented here in a 
number of ways however. For example general symmetric matrices have | n(n — 1) 
degrees of freedom wheras the Kuramoto graphs have n — 1 degrees of freedom (without 
loss of generality any one node can be fixed at phase angle 0). Also we seek negative 
definiteness in conjunction with n equality conditions in the form of the frequency fixed 
point equations (8). 

2.1. The weighted Kuramoto model and the Partition problem 

Using the standard formalism of [10] 

Partition 

Instance: A positive integer m and a collection of positive integers aq, 02 ,.., a m . 
Question: Is there a subset S of {1,2, ..,m} in which 



(13) 


For any collection of m positive integers aq, a 2 , a m we construct a particular weighted 
graph G[ai,a 2 , ..,a m ], shown in figure 2, with corresponding adjacency matrix A t j and 
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Figure 2. The graph G[a 1 , 02 , a m ] 


show that Partition is satisfied if and only if there is a non-zero stable fixed point for 
the corresponding homogeneous Kuramoto model. 

Let t = 1/2 JT a*. Set the vertex set of G to be V = x, u\, 112 ,u m , iq, V 2 , y. 
For each i = 1 , ,.,m there is a link of weight a* between x and Ui, a link of weight a* 
between iq and Vi, and a link of weight 2fa* between v t and y. Let {9*,i = l,..,n} 
be any non-zero stable hxed point of the homogeneous Kuramoto model corresponding 
to Aij. We shall show that there is a non-zero stable hxed point for A t] if and only if 
Partition is satisfied for inputs m, cq, < 22 , a m . To do this we hrst assume that there is 
a subset S satisfying equation (11). Consider the assignment of phase angles shown in 
the phase angle diagram figure 3 where the nodes are assigned from six different phase 
angles shown within the diagram at 0, a, 2a, n, 2n — 2a, and 2n — a. The open nodes 
refer to groups of nodes with the same phase angle. The total edge weights between node 
clusters are indicated by numbers on the links. It is straightforward to verify that this is 
a hxed point solution where a=arcos(l/2t). Since all the phase angle differences across 
edges are also between (— 7t/2,7t/2) the hxed point is also stable. Now assume that the 
system has a non-zero stable hxed point {9*,i = 1, ..,n} of the corresponding weighted 
homogeneous Kuramoto model. Let the stable hxed point angles corresponding to the 
nodes of G be denoted by 9*[x] — 0 at x, 9*[ui\ at Ui, 9*[vi] at Vi and 9*[y\ at y. If for 
each i we apply equation (8) to vertex v t we have 


ta i sm(9*[v i ] - 9*[y\) = a* sin( 6 >*[uj] - 9*[vi\) 


(14) 


so that | sin( 6 , *[n i ] — 9*[y])\ < 1/2 1. For t large this means that 9* [w/ — 9*[y] is close to 
0 or 7T. However by the stability inequality (10) applied at X = {n*}, 


tdi cos^* [vi] — 9*[y ]) + a* 008 ( 6 ** [tij] — 9*[vi\) > 0. 


(15) 


This rules out the latter possibility (that 9*[vi) — 9*[y] is close to 7 r) and for l/2t small 
certainly |6**[uj] — 9*[y\\ < 1 jt. We must now have the phase angles arranged as in the 
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phase angle diagram of figure 4 where / and g are integers with f + g — ^2 a i and say 
/ < g. For stability f3 < 7t/2 so that a > 7t/2 — 1/t. By equation (8) at x 

f >g^{\-\)>g~\- (16) 

Since t > g and / and g are integers we must have f = g and thus we have found 
a partition of ai, a 2 ,a m . Thus we have shown that Partition is satisfied for inputs 
m, ai, a 2 , .., a m if and only if there is a non-zero stable fixed point for the corresponding 
weighted homogeneous Kuramoto model. Since Partition is NP-complete [10] it follows 
that determining whether weighted Kuramoto models have a non-zero stable fixed point 
is NP-hard. Note that it does not follow that this problem is NP-complete since it 
is not at all clear that this problem is in NP. Establishing NP-hardness for the case 
where A is a {0,1} matrix seems to be more difficult. By replacing each vertex u m 
and v m with a clique of size a m we can obtain an unweighted graph G'[ai, a 2 , a m ] 
corresponding to a 2 , a m ] with G' having a non-zero stable fixed point if and only 




if G does. However G' and the corresponding matrix Ah require a i bits of data in 
their representation, while the Partition problem is known to have an algorithm that 
runs in time polynomial in Jj] cq [10]. What we seek is an unweighted graph G such that 
the problem of detecting non-zero stable fixed points might require exponential time in 
the size of G (or is NP-hard in a strong sense [10]). To do this in the following we refer 
to a real valued variation of the Partition problem. 


2.2. The 0-1 Kuramoto model and a partition problem over the reals 


The following describes a variant of the Partition problem over the real numbers which 
we call Kuramoto Partition. In the body of this section we shall construct an unweighted 
graph and corresponding (0,1) adjacency matrix such that Kuramoto Partition is 
satisfied if and only if there is a non-zero stable fixed point for the corresponding 
unweighted homogeneous Kuramoto model. 

Kuramoto Partition 

Instance: A positive integer n and collection of positive integers bi,i = 1, ..,n such that 
2 < b\ < &2 < .. < b n < 7i. 

Question: Let b t = n,i = n + 1,.., 3n. Is there a subset S of {1, 2,.., 3n} in which 


V / *1(1 - W y- / 6?(l~f 2 ) 

k V 1 + b > + 2b “ k V 1 +^ 

1 > e > 0? 

n 

We note by continuity that (17) and (18) have a solution if and only if 

y -A= > y -A= and 

ks V 1 +k a+m 

y I kit 1 ~ 1/"T / ^(1 - l/n 2 ~ 

k V 1 + 6 i + 2b ‘/ n k V 1+6 * ? - 2 b ‘i n 


(17) 

(18) 


(19) 

( 20 ) 


Notice that Kuramoto Partition has input size nlog^n. We show that any polynomial 
algorithm for detecting non-zero stable fixed points in the (unweighted) homogeneous 
Kuramoto model must produce a polynomial algorithm for Kuramoto Partition. We 
briefly provide some background to this type of problem. The Partition problem over 
sets of general real numbers have good heuristic methods. For example the differencing 
method [14] tries to find a partition of a set of real numbers between 0 and 1 that 
minimises the difference between the sum of terms in each part, and runs in time O(n). 
Probabilistic analysis [15] shows that for independent random selections of the real 
numbers, the minimum difference has a median that is exponentially small, while the 
differencing method will find differences as small as 0(n~ cl ° 9n ), c > 0 with the exception 
of pathological or rare examples. This suggests that for most instances Kuramoto 
Partition is satisfied and that solutions can be found in polynomial time. On the other 
hand for some instances Kuramoto Partition is not satisfied (for example b{ = n for 
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Figure 5. The graph G[b 1: b 2 , b 3n } 


all i and n is odd), and is less likely to be satisfied if the bi have few different values. 
In summary despite good heuristic methods the author is not aware of any polynomial 
algorithm (in n ) for this problem nor whether Kuramoto Partition is NP-hard. 

Now let bi,i = 1, ,.,n be integers such that 2 < b\ < 62 < ■■ < b n < n. Also let 
bi = n,i = n + l,..,3n. Construct a graph G[bi, b 2 ,.., b 3n ] as follows. A vertex x is 
adjacent to nodes u±, ..,u 3n , and a vertex y is adjacent to nodes v 3 , ..,v 3n . Additionally 
for each i = 1 ,.., 3n both Ui and Vi are adjacent to separate complete graphs on bi nodes. 
This is illustrated in figure 5. Let A hJ be the adjacency matrix of G and consider the 
corresponding homogeneous Kuramoto model. We shall show that this system has a 
non-zero stable fixed point if and only if the corresponding Kuramoto Partition problem 
is satisfied. To begin we shall assume that we have a solution S to Kuramoto Partition. 
Let the angles 7 , a*, i = 1,.., 3 n be defined by 


7T > 7 = 2arcos(e) 


OLi 


= arsm 


b?( 1~6 2 ) 

1 + fi 2 + 2 bi€ 


= arsm 


7 


6?(1 


1 + bf - 2 b^ 


Pi =^-a.i,ieS 
7 

= 7T - - - Oti,i e S c 



(21) 

i e S 

(22) 

i G S c 

(23) 


(24) 


(25) 


Consider the assignment of phase angles as depicted in the phase angle diagram 
figure 6 where for % e S U{ is located at a t , is located at a* + 2 Pi and all the nodes 
of the adjacent clique are co-located at a* + Pi, where sin a t = bi sin Pi. For * e S c Ui is 
located at 2 tt — a i: v t is located at 2 i r — — 2 Pi and all the nodes of the adjacent clique 
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Figure 6. The phase angle structure of any non-zero fixed point solution 


are co-located at 2tt — oti — /3j. It is a simple matter to verify that 


sinl 


= 

bi sin 


(26) 

7 


= 

2 oti + 2/3for i 

e S 

(27) 

sin( 

Oti) 

= 

h sin (fa), 


(28) 

2 vr- 

~ 7 

= 

2 oti + 2/3,; for i 

e S c 

(29) 

E 

sin( 

[oti) = 

E sin ( a *) 


(30) 

i£S 



ieS c 



cos( 

'2 J 

< 

1 

n 


(31) 

7 


< 

7T 


(32) 


where (30) follows from (17) and (31) follows from (18). Equations (26)-(30) ensure 
that this system is a fixed point solution and (31) ensures that each cq < |. We must 
therefore have a non-zero stable fixed point solution. 

We now assume {0*,i = l,..,n} is any non-zero stable fixed point of the 
corresponding Kuramoto homogeneous model, and show that Kuramoto Partition must 
be satisfied. We do this by showing in stages that the only possible phase angle geometry 
follows the structure of figure 6 and that this in turn leads to a solution to (17)-(18). 
First we show for each k that all of the phase angles of the nodes C k of the clique 
adjacent to u k and v k must all have the same phase angles. Let D k = C k U{ k aJ U{ u fc}- 
Let 


P =j2 z j=Yl e ** = E ( cos ( 0 i)+* sin (0j))- 

j£D k jeD k j&D k 
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Z m Zj 


jm.J 


Then if V is the complete set of all the vertices in G for any m G Ck 
ZmP Zm ^ ^ } ZrnZj 1 T ^ } 

j'eDfe j&D k i£D k ,j^m 

= 1 + COS ( A *jm) +i Sill ( A l 

j&D k ,j^m j&D k ,j^m 

= 1 + 5^ A jm cos(A* m ) + i ^2 A jm sin (A* 
j&v jev 

= l + ]TA jm cos(A*J. (33) 

j£V 

We consider two cases. If p ^ 0 then for some angle 9 P , p = ||p||[cos0 p + isin9 p ). Tims 


jm/ 


ZmP = | |?| | [cOs(9 P '*■ C) + i sin (0 P - 0‘m) ' 


(34) 


Comparing equations (33) and (34) we have sin(d p — 9* J = 0. Also by the stability 
condition (10) cos (9 P — 9* m ) > 0. It follows that 9* n = 9 P . If on the other hand p — 0 
then Y2iev Aim cos(A* m ) = —1 which contradicts the stability condition (10). Thus the 
nodes of every clique in G must also have the same phase angle. In order to relate the 
non-zero stable phase angles to nodes of G we let x be at angle 9*[x] — 0, y at 9*[y\, 
Ui at 9*[uj], Vi at 9*[vi], and the clique of nodes C l adjacent to tq and Vi at 9*[Cj\. For 
each i by the fixed point equation (9) applied with X—Ci we have 

bi sin(0*[uj] - 9*[Ci]) = b t sin( 6 »*[C' i ] - 9*[ui\) = b k d (35) 

for some d between 0 and 1. This means that sin(#*[ui] — 9* [C ,]) and sin(0*[Cj] — 9*[u t ]) 
are either arsin(d) or n — arsin(d). This latter possibility is ruled out however by the 
stability condition (10) applied at X—Ci. Thus 

sin(0*[ui] — 9*[Ci]) = sin( 6 , *[Ci] — 9*[ui\) = arsin(d). (36) 

Similarly by choosing X=Ci UI 1 ’*} then 

sin(^*[y] - 0*[vi\) = sin( 6 , *[rtj] - 9*[x]) < (37) 

The stability condition (10) can also be used to show for each % that the phase angles of 
Vi and C* must all be on the same side of figure 6 in the sense that 9*[ui], 9*[vi\ and 
6* [Q] are between 0 and 9* [y] , or between 6* [y] and 2n. Thus the phase angle geometry 
must follow the structure of figure 6 where we assume by symmetry that 7 < n. The 
fixed point conditions then take the form of equations (26)-(30). Finally we note that 
n + 1, ,.,3n cannot all be in S or all in S c . If for example n + 1, ..,3n G S then for 
j G S c the phase angle differences between y and Vj, and between Uj and x can be 
shown to be greater than | so that the stability inequality ( 10 ) is violated by taking 
A = Ci Ul 'q} Ulw}- The stability condition is therefore equivalent to the condition 
otj < f. By equation (28) and (29) this is satisfied if and only if 7 is at least 7 ' where 
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Thus the stability condition is equivalent to inequality (31). We have now 
established the conditions (26)-(32). By setting e = cos(^) equation (17) follows by 
substitution from (26)-(29) into (30) while (18) follows from (31) and (32). We therefore 
have a solution to Kuramoto Partition. Thus we have shown that Kuramoto Partition 
is satisfied if and only if there is a non-zero stable fixed point for the corresponding 
unweighted homogeneous Kuramoto model. 


2.3. Surd Partition 


The essence of the difficulty of Kuramoto Partition lies in the use of square roots with 
typically infinite decimal expansions. To emphasise this we provide a more generic but 
related partition problem that may be of interest to complexity and number theorists 
and conjecture that this problem is NP-hard. We note that if the square roots are 
omitted in the following then the problem is well known to have a polynomial solution 
[ 10 ]. 

Surd Partition 

Instance: A positive integer n and a collection of positive integers bi,i = 1 ,..,n such 
that 1 < bi < b 2 < .. < b n < n. 

Question: Is there a subset S of {1,2, ..,n} in which 


ieS ies c 


< 1 ? 


(39) 


3. Conclusion 

In this paper we demonstrated the fundamentally complex relationship between the 
network topology and the solution space of the Kuramoto model by showing that de¬ 
termining the possibility of multiple stable fixed points from the network topology is 
NP-hard for the weighted Kuramoto model. Specifically we show that the NP-complete 
Partition problem is satisfied if and only if there is a non-zero stable fixed point solution 
to a related weighted homogeneous Kuramoto model. In the case of the unweighted 
Kuramoto model we show that a particular real partition problem is satisfied if and 
only if there is a non-zero stable fixed point solution to a related unweighted homoge¬ 
neous Kuramoto model. A simplified variant of this partition problem that may be of 
interest to complexity and number theorists is given and conjectured to be NP-hard. As 
a consequence ee conclude that it is unlikely that stable fixed points of the Kuramoto 
model can be characterized in terms of easily computable network invariants. 
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